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Electron transport through a single-level quantum dot weakly coupled to Luttinger liquid leads 
is considered in the master equation approach. It is shown that for a weak or moderately strong 
interaction the differential conductance demonstrates resonant-like behavior as a function of bias 
and gate voltages. The inelastic channels associated with vibron-assisted electron tunnelling can 
even dominate electron transport for a certain region of interaction strength. In the limit of strong 
interaction resonant behavior disappears and the differential conductance scales as a power low on 
temperature (linear regime) or on bias voltage (nonlinear regime). 
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I. INTRODUCTION 

Last years electron transport in molecular transistors 
became a hot topic of experimental and theoretical in- 
vestigations in nanoelectronics (see e.g. 1,2 ). From exper- 
imental point of view it is a real challenge to place a single 
molecule in a gap between electric leads and to repeatedly 
measure electric current as a function of bias and gate 
voltages. Being in a gap the molecule may form chemical 
bonds with one of metallic electrodes and then a consid- 
erable charge transfer from the electrode to the molecule 
takes place. In this case one can consider the trapped 
molecule as a part of metallic electrode and the corre- 
sponding device does not function as a single electron 
transistor (SET). Much more interesting situation is the 
case when the trapped molecule is more or less isolated 
from the leads and preserves its electronic structure. In 
a stable state at zero gate voltage the molecule is electri- 
cally neutral and the chemical potential of the leads lies 
inside the gap between HOMO (highest occupied molec- 
ular orbital) and LUMO (lowest unoccupied molecular 
orbital) states. This structure demonstrates Coulomb 
blockade phenomenon 3,4 and Coulomb blockade oscilla- 
tions of conductance as a function of gate voltage (see 
review papers in 5 and references therein) . In other words 
a molecule trapped in a potential well between the leads 
behaves as a quantum dot and the corresponding device 
exhibits the properties of SET. The new features in a 
charge transport through molecular transistors as com- 
pared to the well-studied semiconducting SET appear 
due to "movable" character of the molecule trapped in 
potential well (the middle electrode of the molecular tran- 
sistor). Two qualitatively new effects were predicted for 
molecular transistors: (i) vibron-assisted electron tun- 
nelling (see e.gAi) and, (ii) electron shuttling 8 (see also 
review 9 - ). 

Vibron(phonon)-assisted electron tunnelling is induced 
by the interaction of charge density on the dot with local 
phonon modes (vibrons) which describe low-energy exci- 
tations of the molecule in a potential well. This inter- 



action leads to satellite peaks (side bands) and unusual 
temperature dependence of peak conductance in resonant 
electron tunnellingiS. For strong electron- vibron interac- 
tion the exponential narrowing of level width and as a 
result strong suppression of electron transport (polaronic 
blockade) was predicte d 10 ! 11 . The effect of electron shut- 
tling appears at finite bias voltages when additionally to 
electron-vibron interaction one takes into account coor- 
dinate dependence of electron tunnelling amplitude 8,9 . 

Recent years carbon nanotubes are considered as the 
most promising candidates for basic element of future 
nanoelectronics. Both Ceo-based and carbon nanotube- 
based molecular transistors were already realized in 
experimen t 12 i 13 . The low-energy features of I-V charac- 
teristics measured in experiment with Ceo-based molec- 
ular transistor 12 can be theoretically explained by the 
effects of vibron-assisted tunnelling 7 . 

It is well known that in single-wall carbon nanotubes 
(SWNT) electron-electron interaction is strong and the 
electron transport in SWNT quantum wires is described 
by Luttinger liquid theory. Resonant electron tunnelling 
through a quantum dot weakly coupled to Luttinger liq- 
uid leads for the first time was studied in Ref. 14 were a 
new temperature scaling of maximum conductance was 
predicted: G(T) oc T x l g ~ 2 with interaction dependent 
exponent (g is the Luttinger liquid correlation parame- 
ter). 

In this paper we generalize the results of Refs j 1Q i 14 to 
the case when a quantum dot with vibrational degrees 
of freedom is coupled to Luttinger liquid quantum wires. 
The experimental realization of our model system could 
be, for instance, Ceo-based molecular transistors with 
SWNT quantum wires. 

In our model electron-electron and electron-phonon in- 
teractions can be of arbitrary strength while electron tun- 
nelling amplitudes are assumed to be small (that is the 
vibrating quantum dot is weakly coupled to quantum 
wires). We will use master equation approach to eval- 
uate the average current and noise power. For noninter- 
acting electrons this approximation is valid for temper- 



2 



atures T ^> r , where r is the bare level width. For 
interacting electrons the validity of this approach (per- 
turbation theory on To) for high-T regime of electron 
transport was proved for g < 1/2 (strong interaction)^ 5 - 
and when 1 — g <C 1 (weak interaction)^. 

We found that at low temperatures: To <C T <C hwo 
(hwo is the characteristic energy of vibrons) the peak 
conductance scales with temperature accordingly to Fu- 
rusaki prediction^: G(T) cx (Y\/T)(T / K) 1 ^ 9-1 (A ~ e F 
is the Luttinger liquid cutoff energy). The influence of 
electron-phonon interaction in low-T region results in 
renormalization of bare level width: Y\ = roexp(— A 2 ), 
where A is the dimensionless constant of electron-phonon 
interaction. In the intermediate temperature region: 
hwo < T < \ 2 hwo, (A ^> 1), Furusaki scaling is changed 
to G(T) cx (T) 1 /^- 3 / 2 and at high temperatures when 
all inelastic channels for electron tunnelling are open we 
again recovered Furusaki scaling with nonrenormalized 
level width (To). 

For nonlinear regime of electron tunnelling we showed 
that zero-bias peak in differential conductance, present- 
ing elastic tunnelling, is suppressed by Coulomb corre- 
lations in the leads. This is manifestation of the Kane- 
Fisher effec t 14 1 15 . When interaction is moderately strong 
(1/2 < g < 1) the dependence of differential conductance 
on bias voltage is non-monotonous due to the presence 
of satellite peaks. For g > 1/2 the zero-bias peak can 
be even more suppressed than the satellite peaks, which 
dominate in this case. This is the manifestation of the 
interplay between the Luttinger liquid effects in the leads 
and the electron-phonon coupling in the dot . For strong 
interaction g < 1/2 satellites are also suppressed and the 
differential conductance at low temperatures (T <C hwo) 
scales as dl/dV cx V x l 9 ~ 2 . This scaling coincides with 
the Furusaki prediction, where temperature is replaced 
by the driving voltage (eV) which becomes the relevant 
energy scale for eV ^> T,hwo,T. It means that the in- 
fluence of vibrons on the resonant electron tunnelling 
through a vibrating quantum dot can be observed only 
for weak or medium strong interaction (1/2 < g < 1) in 
the leads. 



II. THE MODEL 

The Hamiltonian of our system (vibrating quantum 
dot weakly coupled to Luttinger liquid leads, (see Fig.l) 
consists of three parts 



(see e.giii) 



n = nLL + n QD + n T 



(i) 



Here U LL = V . 



=l,r 



describes quantum wires adi- 



abatically connected to electron reservoirs. Quantum 
wires (left-L and right-R) are supposed equal and mod- 
elled by Luttinger liquid Hamiltonians with equal Lut- 
tinger liquid parameters l/g L (R)'- V#L = 1/9r = 1/g 



n 



L(R) 



Hi = hv c / atakkdk . 



(2) 



Here a^(ak) are the creation (annihilation) operators of 
bosons which describe the charge density fluctuations 
propagating in the leads with velocity v c ~ vp. These op- 
erators satisfy canonical bosonic commutation relations 
[afc,ajj",] = S(k — k'). In what follows we consider for 
simplicity the case of spinless electrons. 

The Hamiltonian of vibrating single level quantum dot 
takes the form (see e.g. 10 ) 

n QD = eof+f + e % (b + + &)/+/ + hw b+b , (3) 

where Sq is the energy of electron level on the dot, Hwq is 
the energy of vibrons, Si is the electron- vibron interaction 
energy, f + (f) and b + (b) are fermionic (/) and bosonic (b) 
creation (annihilation) operators with canonical commu- 
tation relations {/, / + } = 1, [6, b + ] = 1. 

The tunnelling Hamiltonian is given by standard ex- 
pression 



(4) 



j=L,R 



where tj is the electron tunnelling amplitude and ^f(j), 
j = L, R is the annihilation operator of electron at the 
end point of L(R)-electrode. This operator could be writ- 
ten in a "bosonised" form (according toi^) 



*(£(#)) = \ — -exp 



Jo 



dk- 



- oik 1 2 



4) 



(5) 

Here a is a short-distance cutoff of the order of the recip- 
rocal of the Fermi wave number kp and K p = (2/ g — l)~ l 
is the interaction parameter in the "fermionic" form of 
the Luttinger liquid Hamiltonian (2), it defines the Lut- 
tinger liquid parameter g which is varied between and 
1: the case g = 1 describes the " noninteracting" (Fermi- 
liquid) leads, than in the case g — » 1 the interaction in 
the leads goes to infinity. 

Hamiltonian (3) is " diagonalized" to 
H>d = £pf + f + hwob+b by the unitary trans- 
formation (see e.gJ^) U = exp(zApnj), where 
p — i(b + — b)/y/2, rif = f + f and the dimensionless 
parameter A = —\/2Ei/hwo characterizes electron- vibron 
coupling. The unitary transformation results in: (i) the 
shift of fermionic level (polaronic shift) Sp = Sq — e 2 /hwo 
and, (ii) the replacement of tunnelling amplitude in 
(3) tj => tj • exp(-iA • p). The model Eqs.(pQ)-([5j)can 
not be solved exactly and one needs to exploit certain 
approximations to go further. 

We will use "master equation" approximation (see 
e.g. 5 ) to evaluate the average current and noise power in 
our model. It is in this approximation that average cur- 
rent separately for the model with interacting leads 14 and 
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FIG. 1: The schematic picture of the two-terminal electron 
transport through a vibrating quantum dot weakly coupled 
(via narrow dielectric regions Ht) to the Luttinger liquid leads 
(Tif^ — Hi) with the chemical potentials £F±eV/2 (V is the 
driving voltage) . All the energies are counted from the Fermi 
energy, which chosen to be zero. Electrons tunnel from one 
lead to another by hopping on and off the dot level with the 
energy ep (elastic channel) and due to electron- vibron cou- 
pling they can emit or absorb vibrons (vibron-assisted tun- 
nelling). Inelastic channels are represented as side- levels with 
energies ep± =£p± huuo. The position of the dot levels with 
respect to the Fermi energy can be uniformly shifted by ap- 
plying voltage V g to the "gate" electrode. 



for vibrating quantum dot with noninteracting leads^ 
was calculated earlier. Master (rate) equation approach 
exploits such quantities as the probability for electron 
to occupy dot level and the transition rates. It neglects 
quantum interference in electron tunnelling and there- 
fore describes only the regime of sequential electron tun- 
nelling which is valid when the width of electron level 
To <C min(T, eV) . In other words, in our case "master 
equation" approach is equivalent to the lowest order of 
perturbation theory in Tq. 

For interacting electrons the validity of master equa- 
tion approach for high-T regime of resonant electron 
tunnelling can be justified for strong repulsive interac- 
tion g < It is correct also for weak interaction 
1 — g <C 1 as one can check by comparing the results of 
Refjii and Ref^, where resonant tunnelling through a 
double-barrier Luttinger liquid was considered for weak 
electron-electron interaction. Notice, that the results 18 
of exact solution known for g = 1/2, where a mapping 
to free-fermion theory can be used 5 , do not agree with 
the high-T scaling of G(T)^ extrapolated to this special 
point g = 1/2. The free-fermion scaling G(T) oc T _1 
found for g = 1/2 (master equation approach predicts 



T-independent valued) could be a special feature of this 
exactly solvable case. We will assume that beyond the 
close vicinity to g = 1/2 the master equation approach 
for high-T behavior of conductance is a reasonable ap- 
proximation. 



III. TRANSITION RATES AND THE AVERAGE 
CURRENT 

In master equation approach the average current 
through a single level quantum dot expressed in terms 
of transition rates takes the form 



~pR rL pL p 

, i 01 i 10 1 01 1 



10 



(6) 



where is the rate of electron tunnelling from the 

R( L) 

dot to right (left) electrode, r l0 describes the reverse 
process and T s = r i +r 10 , T if = rf f +T? f (i, f = 0, 1). 
To evaluate these rates in our approach we will use Fermi 
"Golden Rule" (quantum mechanical perturbation the- 
ory) for tunnelling Hamiltonian obtained from Eq.(4) af- 
ter the unitary transformation: Ht => Ht 

Kt = J2 {*i* + 0")/expHA.p) + ft.c.} . (7) 

j=L,R 

The standard calculation procedure results in the follow- 
ing expressions for tunnelling rates 
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1 oi — 



/oo 
dt(V(t)V + (0)) b (*+(t)*j(0)) f 
-oo 

• exp(z(ep — sf + eVj)t/h) , (8) 



dt(v + (t)v(o)) b (* 3 mi(o)) f 

• exp(— i(sp — sf + eVj)t/h) , (9) 



where Vl — Vr = V is the bias voltage and j = L,R. 
Notice that in the perturbation calculation on the bare 
level width To ex l^,^! 2 , we neglect the level width in the 
Green function of the dot level. Besides, in this approxi- 
mation averages over bosonic and fermionic operators in 
formulas for tunnelling rates are factorized and, thus, the 
averages (. . over bosonic variables 



V = exp(— iX • p) 



b] 



(10) 



can be calculated with the quadratic Hamiltonian 
Hb = hwob+b. In what follows we will assume that vi- 
brons are characterized by equilibrium distribution func- 
tion: rib{T) = [ex.-p(hwo/T) — Averages (...)/ over 
fermionic operators in Eqs.(8),(9) are calculated with the 
Luttinger liquid Hamiltonian (2). The corresponding cor- 
relation functions in Eqs.(8),(9) are well known in the 
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literature (see e.g ji 14 ) 



(V(t)V+(0)) b = exp(-A 2 (l + 2n B )) 



K I— — 00 



■I J ' [2AV«b(1 + n B )\ exp[-ilw (t + iK/2T)\ 



(11) 



voltage and its value is totally determined by the effective 
level width Y. 

For a vibrating quantum dot (A ^ 0) weakly coupled 
to noninteracting leads (g = 1) our approach reproduces 
the results of Ref. 10 . In the temperature region we are 
interesting in (T ^> r^m) the general formulae derived 
in 10 can be rewritten in a more clear and compact form. 
In particular, by using for g = 1 in Eq.(ll) the well- 
known representation for Gamma function (see e.g. 19 ) 



A 



iA_ 
irT 



sinh 



ttT -t 



h 



-1/9 



(12) 

Here Ii(z) is a modified Bessel function, A ~ £p is a ultra- 
violet cutoff energy, g is the Luttinger liquid correlation 
parameter. 

By putting correlation functions (11), (12) in 
Eqs.(8),(9) and evaluating time integrals we get 

<j) 

10 



the following equations for tunnelling rates 



and T 



r%,(j = L,R) 



1 a exp [-A 2 coth(fi^ /2T) ± Aj/2T] 



2irh 
/2ttT 



V A 



1/g-l oo 



A 2 



_smh(hw /2T)_ 



(13) 



where r^^) is the partial level width (see, for 
example,— ) Tj = {^ncjt 2 - / ah) {f + f) = const j 
(j = L,R), T L + T R = T , A, = e F - e P + eV-; here 
r(z) is Gamma function. 

At first we consider different limiting cases when it 
is possible to obtain simple analytical expressions for the 
average current (6). Notice, that electric current depends 
on the gate voltage V g through the corresponding depen- 
dence of level energy Sp(V g ). It is convenient for the fur- 
ther analysis to choose the value of gate voltage at which 
the current at low bias is maximum as: Sp(V g ) = Sf- In 
what follows we also put Vl = —Vr = V/2. 

For noninteracting leads (g = 1) and noninteracting 
quantum dot (A = 1) it is easy to derive from Eqs.(6),(13) 
the well-known formula for the maximum (resonant) cur- 
rent at temperatures T ^> ^l(r) 



t /-t t\ eT 1 [eV\ 
W- Y tanhf-j , 



(14) 



where V = T a/ (Y F + Tr) is the effective level width. 
It is evident that at high voltages: eV ^> T the current 
through a single level dot does not depend on the bias 



• iz 



COSh(7Tz) 



(15) 



it is easy to obtain the following expression for the max- 
imum (peak) conductance 



G X (T) = G{T)-F X 



( hw \ 



where 



7T r e 2 
G(T)c±-G - , G = - 



(16) 



(17) 



is the standard resonance conductance of a single-level 
quantum dot at T ^ ^l(r)- The dimensionless function 
F\(x) takes the form 



F x (x) =exp[-A 2 (l + n B (x))] 

h{z{x)) \ 
^ cosh 2 (lx/2) J ' 



(18) 



here z(x) = 2\ 2 ^riB{x)[l + tib(x)] and 
tib(x) = (exp(x) — l) -1 . At low temperature re- 
gion r <C T <C hwo, when there are no thermally 
activated vibrons in the dot (n# <C 1) only the first 
term in the brackets contribute to the sum and: 
F\(T <C hwo) ~ exp (—A 2 ). We see, that zero-point 
fluctuations of the dot position result in renormal- 
ization of the level width T\ = T • exp (—A 2 ). For 
strong electron-vibron coupling this phenomenon (po- 
laronic narrowing of level width) leads to polaronic 
(Franck- Condon) blockade of electron transport through 
vibrating quantum dot^. The temperature behavior of 
peak conductance (16) was considered in Ref^. 

Now we will study the general case when interact- 
ing quantum dot (A ^ 0) is connected to interact- 
ing leads (g < 1). Analytical expressions for conduc- 
tance in this case can be obtained in the limits of low 
(^L(R) «^« hwo) and high (T ^> Hwq) temperatures. 

At low temperatures the main contribution to the sum 
over "Z" in Eq.(13) comes from elastic transition I = 0. 
All inelastic channels (I ^ 0) are exponentially sup- 
pressed for eV, T <C Hwo. At T <C ghwo the peak con- 
ductance takes the form 



G(T) * ^Go^ 



r(i/2 5 ) 



|T (1/2$ + 1/2) 



TTT' 

A 



i/s-i 



(19) 
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We see that at low temperatures conductance scales 
with temperature according to Furusaki's prediction^: 
G(T) rsj T 1 /^ -2 . The influence of electron- vibron cou- 
pling results in multiplicative renormalization of bare 
level width F\ = Texp(— A 2 ). 

At high temperatures: T ^> Hwo one can use the 
well known asymptotic expansion for Bessel function 
Ii(z) c± exp(z)/V27rz, which can be used in summation 
Eqs. ([T8]h (13) until I 2 < z. Besides, in this temperature 
region the summation in Eq.(13), can be replaced by in- 
tegration and the corresponding integral can be taken 
exactly 



f 

j — c 



\T(a + iz)\ 2 dz = 7r2 1 " 2a r(2a) . 



(20) 



This allows us to derive the following expression for 
the temperature dependence of peak conductance in the 
intermediate temperature region hwo <C T < \ 2 hwo, 
(A>1) 



G(T) ~ -Go 



rexp {-\ 2 hw /4T) 



(21) 

Notice that in the considered temperature region 
the polaronic blockade is already partially lifted 
r A (T) = Texp (-\ 2 hw /4T) ~ T at T ~ X 2 hw and con- 
ductance scales with temperature as G(T) ~ T x l 9 ~ z l 2 . 
At last, at temperatures T ^> X 2 hwo when all inelastic 
channels for electron transport are open, the polaronic 
blockade is totally lifted 20 and we reproduce again Fu- 
rusaki scaling. It is clear from our asymptotic formulae 
(19), (21) that both in low- and in high- temperature re- 
gions the contributions of electron-electron and electron- 
vibron interactions to the conductance are factorized. In 
general case these contributions are not factorized, as one 
can see from Eqs.(8),(9) and from Eq.(13) for tunnelling 
rates, and we can expect interplay of Kane-Fisher effect 
and the effect of phonon(vibron)-assisted tunnelling. 

To see this interplay we consider nonlinear (differen- 
tial) conductance G(V) = dl/dV. It is well known 
that Kane-Fisher effect is pronounced for the energies 
close to the Fermi energy. For differential conductance 
it means that the zero-bias (elastic) resonance peak is 
suppressed with the increase of electron-electron inter- 
action, while satellite peaks are less affected by the in- 
teraction. When electron-electron interaction is weak or 
moderately strong (1/2 < g < 1) the dependence of dif- 
ferential conductance on the bias voltage (for A ~ 1) is 
not a monotonous function due to the presence of satel- 
lite peaks (see Figs. 2, 3). 

The resonant behavior disappears for strong inter- 
action g < 1/2 (Figs. 2, 3), when at low temperatures 
T <C hwo differential conductance scales with bias volt- 
age as Giy) oc y 1 /^ -2 in accordance with the Lut- 
tinger liquid prediction for nonlinear electron transport 
through a single- level quantum dot. For instance, if we 
put 1/g = n, (n = 2, 3, . . .) and tune the level energy 




FIG. 2: Differential conductance (in the units of Go) as a 
function of driving voltage (in the units of hwo) for the case 
g > 1/2. Here we put F/kT = 0.01; hw /kT = 10; A 2 = 1; 
and tune the level energy to the resonant position £q — X 2 hwo 
(ep = 0). Solid line corresponds to the case of noninteracting 
leads: g = 1. Dot line corresponds to value g = 0.8, while 
dash-dot line - to g = 0.6. Figure shows how zero-bias (elas- 
tic) resonance peak is gradually suppressed with the increase 
of electron-electron correlations (decrease of Luttinger liquid 
parameter g) while the satellite peaks survive until g > 1/2. 




FIG. 3: Differential conductance (in the units of Go) as a 
function of driving voltage (in the units of hwo) for the case 
g < 1/2. Here F/kT = 0.01; hw /kT; A 2 ; and £o are the same 
as on Fig. 2. Solid line corresponds to the case of noninter- 
acting leads: g = 1; dot line corresponds to value g = 0.45; 
dash-dot line - to g = 0.25. One can see from the figure that 
for g < 1/2 the resonance- like behavior of differential con- 
ductance disappears and conductance scales as a power-law 
of the bias voltage. 



to the resonance point Sq = \ 2 Hwq -("resonant" location 
of the level in the presence of "polaronic" shift), we ob- 
tain for the differential conductance G(V) the following 



FIG. 4: Differential conductance (in the units of Go) as a 
function of level energy eo, counted from the Fermi energy for 
the case g > 1/2. Solid line corresponds to the case of nonin- 
teracting leads: g = 1; dot line corresponds to value g = 0.8; 
dash-dot line - to g = 0.6. The bias voltage eV/ftwo = 5 
is sufficiently high to excite vibrons and to support electron 
transport through inelastic channels. All parameters are the 
same as for Fig. 2. 



FIG. 5: Differential conductance (in the units of Go) as a 
function of level energy eo, counted from the Fermi energy for 
the case g < 1/2. Solid line corresponds to the case of nonin- 
teracting leads: g = 1; dot line corresponds to value g = 0.45; 
dash-dot line - to g = 0.25. The bias voltage eV/ftwo = 5 
is sufficiently high to excite vibrons and to support electron 
transport through inelastic channels. All parameters are the 
same as for Fig. 3. 



expression for eV/Hwo ^> (1/g — 2) 

1 



G(V)~AkgJ^- 



A l(n-l)\_ 



eV 
2A 



n-2 



(22) 



where T\ = Texp(— A 2 ). One can readily see that ex- 
pression (22) reproduces Furusaki temperature scaling 
Eq.(19) when eV is replaced by T. 

Analogous interplay of Kane-Fisher and polaronic ef- 
fects one can see on Figs. 4, 5, where differential conduc- 
tance is plotted as a function of level energy eo (or, equiv- 
alently, as a function of gate voltage). For noninteracting 
leads (g = 1) the resonance conductance peaks corre- 
spond to the level positions at Ep± = Ef ± eV/2 (in our 
plot we put: eV = 5Hwq). This elastic resonance peak is 
suppressed by electron-electron interaction in the leads 
(g < 1). The dependence G(V g ) for weak and moder- 
ately strong interaction still reveals resonance structure 
with 4 satellites in our case (see Fig. 4). The inelastic 
resonance peaks disappear at g < 1/2 and maximum of 
differential conductance corresponds at g <C 1 to the level 
position at £p(V g ) = £p, that is exactly in the middle be- 
tween chemical potentials of left and right electrodes (see 
Fig.5). 

It is important to stress here once more that for mod- 
erately strong electron-electron interaction in the leads 
the inelastic tunnelling can dominate in electron trans- 
port. One can see from Figs. 2-5 that there is region of 
coupling constants when the first satellite peak is higher 
than the "main" (zero-bias) resonant peak, which corre- 
sponds to elastic (I = 0) tunnelling channel. It is the 
most significant prediction, we have made in this paper. 



IV. THE NOISE POWER 

The knowledge of tunnelling rates Eqs.(8,9) allows us 
to evaluate not only the average current Eq.(6) but the 
noise power as well. We will follow the method developed 
in Refs i 21 i 22 where quantum noise was calculated for res- 
onant electron transport through a quantum dot weakly 
coupled to noninteracting electrodes. 

The noise power is defined (see e.g. 23 ) as 



/oo 
dtexp(iwt)(AI(t)AI(0)) , 
-OO 



(23) 



where AJ(t) = I(t) — I (I is the average current). The 
noise defined in Eq.(23), in the case of sequential tun- 
nelling through a quantum dot, can be expressed in terms 
of tunnelling rates. For a single level quantum dot this 
formula for low frequency noise S = S(w = 0) takes the 
form 

4J 2 



S = 2el 



4e 



2 L 01 1 10 



(24) 



here the average current I is determined by Eq.(6). 
The noise power Eq.(24) depends on temperature and 
bias voltage S(T, V) and contains both thermal (Jonson- 
Nyquist) noise S JN (T) = S(T,V = 0) = 4TG(T) (G is 
the conductance) and the non-equilibrium (shot) noise 
S 3 h(T,V). Since the thermal noise is totally deter- 
mined by temperature dependence of conductance, we 
will study in what follows only shot noise and Fano fac- 
tor F = S s h/2el. In particular, Fano factor in our case 
can be represented as follows 



21 



2e 

T 



•nL pi? 
1 01 1 10 



(25) 
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For noninter acting leads (g = 1) and noninter acting 
quantum dot (A = 0) one readily gets from Eqs.(13),(24) 
a simple expression for the "full" noise (S) of a single 
electron transistor (SET). On resonance £p(V g ) = £f 
and at temperatures T ^> T one finds 



2e 2 r feV 
— - — tann — — 
h \4T 



1 - 

h 



2T feV 
— - tann — — 

'exp(-eF/4T) 



cosh^(eF/4T)_ 



(26) 



From Eq.(26) in the limit V —> we obtain 
S = S JN (T), where Sjn(T) = e 2 T/h = 4TG is the 
thermal noise. In the opposite case eV > T we rederive 
the well-known formulae for the shot-noise and the Fano 
factor of a single level quantum dot* 2 ^ - — 



Ss 



2e 2 r 

h 



2T 



F = 1 



2T 



(27) 



These formulae (26), (27) can be also re-derived from 
the general expression for the full noise of noninteracting 
electrons (see e.g., Eq.(61) in Ref. 23 ) 

S(V,T) = ^J deT t {e)[f L {\ - f L ) + f R (l - f R )} 

+ e l J deTt {e) [1 - T t {e)\ (f L - f R f , (28) 

where T t (e) is the transmission coefficient and 
fj = {exp [(e — Hj)/T] + 1}~ is the equilibrium distri- 
bution function of electrons in the leads (/ij is the chem- 
ical potential; j = L,R). In the case of single level quan- 
tum dot T t (e) takes the form Breit-Wigner tunnelling 
probability 



Tt(e) 



Li- R 



{e-e P Y + {T L +T R f/A 



(29) 



where V = TlTr/ (Tl + I\r). For a weak tunnelling 
when ^l(r) are the smallest energy scales in the prob- 
lem the Lorentzian shape of the Breit-Wigner resonance 
shrinks to (5-function 



T t (e)\r L , R ^ ^2irT6(e-ep) . 



(30) 



With the help of Eqs.(28),(30) for the resonance con- 
dition ep(V g ) = £f we easily re-derive Eqs.(26). (No- 
tice, that in sequential tunnelling approach the tunnelling 
transitions through the left and right barriers are as- 
sumed to be weak and uncorrelated. Therefore we can 
safely neglect T 2 -term in Eq.(28)). It is evident from 
Eqs.(25),(27) that for noninteracting electrons the Fano 
factor is sub-Poissonian (F < 1) and F approaches 1 
for strongly asymmetric junction ^l(r) ^> ^R(L) an d for 
cV > T. 

The master equation approach we have used in our 
analysis holds when electron tunnelling amplitudes are 



1 dS 



sh 



eG dV 




/Tin 



FIG. 6: Differential shot noise power (in the units of eGo) 
as function of the level energy eo in the nonlinear transport 
regime (when eV/huoo = 5) for the case g < 1/2. Solid line 
corresponds to the case of noninteracting leads g = 1; dot line 
corresponds to value g = 0.45; dash-dot line - to g = 0.25. 
Other parameters are the same as on Fig. 5. 



F o 




So/flWo 

FIG. 7: Fano factor for the case g < 1/2 as the function of 
the level energy eo in the nonlinear transport regime (when 
eV/ftwo = 5) which corresponds to Figs. 5, 6. Solid line cor- 
responds to the case of noninteracting leads g = 1; dot line 
corresponds to value g = 0.45; dash-dot line - to g = 0.25. 
All other parameters are the same as on Figs. 5, 6. 



small. For noninteracting electrons this assumption is 
satisfied when electron energies are far from the resonant 
energy level, i.e. T t (e) <C 1. The differential shot noise 
in this case as a function of bias voltage or gate voltage 
behaves similarly to the differential conductance. Notice 
however that due to different dependence on temperature 
the shot noise unlike the thermal one even in sequential 
tunnelling regime (T ^> T) can not be expressed in terms 
of conductance. 

By comparing Fig. 6 with Fig. 5, one can see that the 
above similarity is preserved for interacting electrons 
(g 7^ 1, A ^ 0) as well. The corresponding Fano factor 
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which is the "shot noise/current" ratio and thus is less 
sensitive to the details of tunnelling process, for strong 
electron-electron interaction exhibits a simple behavior 
(see Fig. 7). It dips (F ~ 1/2) at symmetric (with re- 
spect to chemical potentials of the leads) position of the 
dot level. Outside this region F —> 1 (Poissonian noise). 
The width of the dip decreases with the increase of in- 
teraction (see Fig. 7). 

V. SUMMARY 

We considered the influence of interaction on transport 
properties of molecular transistor which was modelled as 
a vibrating single-level quantum dot weakly coupled to 
the Luttinger liquid leads. We found interesting inter- 
play between polaronic and Luttinger liquid effects in our 



system. In particular it was shown that for weak or mod- 
erately strong interaction (1/2 < g < 1) the differential 
conductance demonstrates resonance-like behavior and 
for moderately strong interaction inelastic channels can 
even dominate in electron transport through a vibrating 
quantum dot. For strong interaction < 1) the res- 
onant character of vibron-assisted tunnelling disappears 
and the differential conductance scales as a power law on 
temperature (linear regime T ^> eV) or on bias voltage 
(nonlinear regime eV ^> T). 
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